The gluon spin contribution to the proton spin from the double hehcity asymmetry in 
inclusive production in polarized p + p collisions at ^/s = 200 GeV 
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The double helicity asymmetry in neutral pion production for pt =1 to 12 GeV/c has been 
measured with the PHENIX experiment in order to access the gluon spin contribution, AG, to the 
proton spin. Measured asymmetries are consistent with zero, and at a theory scale of /i^ = 4 GeV'^ 
give AG'""^'"-^' = 0.1 to 0.2, with a constraint of -0.7 < AG'° °^'°-^' < 0.5 at Ax^ = 9 (~3a) 
for our sampled gluon momentum fraction {x) range, 0.02 to 0.3. The results are obtained using 
predictions for our measured asymmetries generated from four representative fits to polarized deep 
inelastic scattering data. We also consider the dependence of the AG constraint on the choice of 
theoretical scale, a dominant uncertainty in these predictions. 

PACS numbers: 13.85.Ni,13.88.+c,21.10.Hw,25.40.Ep 



The quark spin contribution to the proton spin was 
found to be only ~25% P, i, i, indicating that the ma- 
jority of the proton spin on average comes from the gluon 
spin contribution, AG, and/or from gluon and quark or- 
bital angular momentum. High energy polarized proton- 
proton collisions at the Relativistic Heavy Ion Collider 
(RHIC) at Brookhaven National Laboratory access AG 
at leading order through spin-dependent gluon-gluon {gg) 
and quark-gluon {qg) scattering. 

This paper presents results from the 2006 RHIC run 
(Run-6) on AG from measurements of the double helicity 
asymmetry (All) in inclusive mid-rapidity tt" produc- 
tion by the PHENIX experiment. AG can be extracted 
from All using next to leading order (NLO) perturba- 
tive quantum chromodynamics (pQCD) 0], which suc- 
cessfully describes unpolarized cross-sections measured 
at RHIC for many inclusive processes [1, 0, Q , including 
mid-rapidity 7r° production [8|, at ^/s=2QQ GeV. These 
data represent a factor of two improvement in the statis- 
tical uncertainty compared to previous results 
They significantly constrain AG, as presented in a re- 
cent global fit (DSSV) [Til of both RHIC and polarized 
deep inelastic scattering (pDIS) data, which used a pre- 
liminary version of these results. We further present the 
impact of experimental systematic and several theoreti- 
cal uncertainties on our determination of AG. 

We define AG'"'''' = jj^ dxAg{x, fi^), with Ag{x,fi'^), 
the polarized gluon distribution, a function of x, the 
gluon momentum fraction and /i^, the factorization scale. 
Thus AG'"'^! = AG. Figure[TJi shows the best fit Ag{x) 
from four different pDIS fits: GRSV-std [13], BB ("ISET- 
4") [H, LSS [li] and GS-C O] which assumes a node, 
or sign change, in Ag{x). As the pDlS data have lim- 
ited sensitivity to AG, there remains large uncertainty, 
which was estimated by BB, and is shown as the hatched 
band. The result from the recent global fit, DSSV, is also 
shown. It has a node, which is driven by the inclusion 
of RHIC PHENIX 7r° and STAR jet [H [13] All data 



along with evolution from pDIS at large x [ll| . Table |T] 
lists AG'"'''' for several x ranges. 

We define A^l = (o'-i-i-^o'^ )/(o-++ + (t^ ), with o-++ 

(cr^ ) the beam helicity dependent differential cross sec- 
tions for inclusive tt*' production from collisions of lon- 
gitudinally polarized protons with the same (opposite) 
helicity. The asymmetry is measured using 



1 



(PbPy) N++ 



R = 



(1) 



where Pb and Py are the polarizations of the two RHIC 
beams, called "Blue" and "Yellow," and R, the relative 
luminosity, is the ratio of integrated luminosities (L) for 
same and opposite helicity collisions. Here we take N to 
be the tt" yield in a transverse momentum (px) bin. 

Each tt" pt bin samples a distribution in gluon x. Fig- 
ure [T]d shows the sampled gluon x distributions for three 
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FIG. 1: (color online) (a) The polarized gluon distribution as 
a function of x for five fits to polarized data. Hatched band 
is pDIS 1(T uncertainty (BB). (b) Distributions of gluon x in 
three tt^ pt bins from a NLO pQCD simulation. 
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FIG. 2; (color online) (a) Asymmetry in tt" production as a function of pt- Error bars are statistical uncertainties. An 8.3% 
scale uncertainty due to the uncertainty in beam polarization is not shown. The pT independent uncertainty of 7 x 10"'' due 
to relative luminosity is shown only in the inset as a shaded bar. For comparison, we also show our Run-5 result. NLO pQCD 
expectations based on different inputs for AG in the GRSV parameterization are plotted, (b) The profile as a function of 
ACqj^Iy "^' using the combined Run5 and Run6 data considering only statistical uncertainty, or also varying by zblcr the two 
primary experimental systematic uncertainties, from beam polarizations (±AP) and relative luminosity (±A_R). 



All PT bins from a NLOpQCD simulation They 
are peaked at xt/0.7 [10|, with xt = pt/(\/s/2). The 
bins overlap, with our data covering primarily the range 
0.02 <x < 0.3, and so we probe AG^°-°^'°-^l 

The highly segmented PHENIX electromagnetic 
calorimeter (EMCal) [lOl is used to detect tt" 77 
decays. The EMCal covers a pseudorapidity range of 
I77I < 0.35 and azimuthal angle range of A(j) = tt, with 
segmentation At] x A(/) = 0.01 x 0.01. We required 
for each of the two decay photons an energy deposition 
pattern consistent with an electromagnetic shower, no 
charged track pointing to the location of the deposited 
energy, and standard quality assurance requirements Q . 
Events were obtained from an EMCal based high px pho- 
ton trigger [2l| in coincidence with a minimum bias trig- 
ger Q (also used to obtain the relative luminosity) . This 
EMCal based trigger had an efhciency for tt" of 5% at 
PT « 1 GeV/c and 90% for > 3.5 GeV/c. The mini- 
mum bias trigger was defined as the coincidence of signals 
from forward and backward beam-beam counters (BBC) 
with full azimuthal coverage located at pseudorapidities 
±(3.0—3.9) The analyzed data sample corresponded 
to an integrated luminosity of 6.5 pb~^. 

Each collider ring of RHIC was filled with up to 111 
out of a possible 120 bunches, spaced 106 ns apart, 
with bunch helicities set such that all four beam he- 
licity combinations occurred in sequences of four bunch 
crossings. The pattern of helicity combinations for each 
RHIC fill (typically 8 hrs) was cycled between four pos- 
sibilities in order to reduce systematic uncertainties that 



could be correlated to the bunch structure in RHIC @]. 
Events were tagged with the bunch crossing number to 
obtain the beam helicities for the event. The luminos- 
ity weighted beam polarization product was (PbPy) = 
0.322 ± 0.027 (8.3%), with single beam polarizations of 
0.560 and 0.575. Using very forward neutron production 
asymmetry [2^ , the longitudinal polarization fractions 
were found to be greater than 99%. 

As in our previous analyses [1, Q , the relative luminos- 
ity ratio R was obtained from crossing-by-crossing col- 
lected minimum bias (BBC) trigger counts, which mea- 
sure about half of the p+p inelastic cross section [21] . The 
uncertainty on R was derived from the comparison with 
a second trigger based on the Zero Degree Calorimeters 
[23 |. which selects different physics processes in a differ- 



TABLE I: AG for different x ranges at /i^ = 4 GeV^ for each 
group's best fit and the when comparing the expected All 
in Fig. |3ja) with our data (8 degrees of freedom). Also, 
the minimum Bind corresponding /\(j[o o2,o.3l fQ^jj^j jjj pjg 
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ent kinematic range. It contributed a pT independent 
systematic uncertainty to ^ll of 7 x 10^''. 

Equation [T] is used to determine, on a fill by fill basis, 
^LL for the yield in the 7r° mass peak for each pT bin. 
The asymmetries were averaged over fills and corrected 
for the asymmetry in the background contribution (de- 
termined from two 50 MeV/c^ wide sidebands on either 
side of the 7r° peak) [9j , which was consistent with zero. 

Figure [2^ shows the measured from Run-6 [2^ in 
comparison with our published data from the 2005 RHIC 
run (Run-5) Q . The results are found to be statistically 
consistent with a 13% confidence level. The inset shows 
an expanded view of the lowpx region, as well as the rela- 
tive luminosity uncertainty from Run-6. Besides this and 
the scale uncertainty from polarization, other systematic 
uncertainties that can be found by using a bunch polar- 
ization sign randomization technique and by varying the 
TT*^ identification criteria 9] appear negligible. 

Also shown in Fig. [2^ are NLO pQCD predictions of 
All II based on fits of pDIS data by GRSV with three 
different values for AG at the input scale of /x^ = 0.4 
GeV^: 1) "std", their best fit value with AG = 0.24, 2) 
AG = and 3) AG = —1.05. The measurements are 
most consistent with GRSV AG = 0. CTEQ6 unpolar- 
ized parton distribution functions (PDF) |26| were used, 
along with DSS fragmentation functions^ (FF) |2' H|. i n all 
calculations. Using alternative PDF |28|] or FF [23] did 
not lead to significant differences in the All expectations. 

All expectations based on fits to the pDIS data with 
a range of inputs for AG[°'^1 evolved to fi^ = 0.4 GcV^ 
in the GRSV parameterization were calculated. Similar 
to our previous analysis @, values were calculated 
using our combined Run-5 and Run-6 data (25| for these 
expectations, effectively fitting AG with our data in this 
parameterization. In Fig. [2)3, these values are plotted as 
a function of AG^°-°^'°-^^ at ^ 4 GeV^ . Due to soft 
physics contamination at low p^ , we use only data with 
Pt > 2 GeV/c [1]. Assuming that jj, = p^, = A GeV^ 
is then the minimum cutoff of our data. The solid curve 
shows the result considering only statistical uncertainties. 
Due to gg interactions \up+p collisions, All probes AG 
quadratically in 7r° production Q . The profile is thus 
not parabolic, and so we show Ax^ = ~ Xmin = 1 ^^"^ 
9 corresponding to "Icr" and "3cr" uncertainties. The gg 
scattering increases toward low px, which dominates our 
statistics and causes the two minima seen in Fig.[2]3. The 
larger allowed negative region arises from cancelation of 
gg and qg terms in A^l when AG is negative. 

For a robust interpretation of our results in terms of 
AG, we consider not only statistical but also experi- 
mental systematic and theoretical uncertainties. The 
effects of the two largest experimental systematic un- 
certainties, due to polarization and relative luminosity, 
are shown in Fig. ^p. The polarization uncertainty is 
insignificant when extracting AG. However, the uncer- 
tainty on relative luminosity, though small, cannot be 
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FIG. 3: (color online) (a) tt" asymmetry expectations for dif- 
ferent l^gix) in Fig.[T]3. Hatched band is the pDIS uncertainty 
(BB). Combined Run-5 and Run-6 results are also plotted 
(statistical errors only), (b) The profile as a function of 

^(.[0.02.0.3] fo^ 

same parameterizations. Arrows indicate Icr 
uncertainty on BB best fit. Ax^ values are shown for GRSV. 



neglected. Accounting for statistical uncertainty, we find 
^G^^sw'^^ = 0.2 ± 0.1 (Icr) and 0.2l°j (3cr) with an 
additional experimental systematic uncertainty of ±0.1. 

Figure [3^ shows All expectations [l,[l3| based on the 
parameterizations discussed above, along with the pDIS 
uncertainty on ^g{x) in BB propagated to All- Our 
combined Run-5 and Run-6 tt'^ results [1^ are also plot- 
ted for comparison. The values for comparing each 
curve with our data are given in Table [T| The three fit re- 
sults without a node in Ag(x)-LSS, GRSV and BB-have 
large values of AG[° °^'°''^1 which lead to relatively large 
asymmetries that lie mostly above the data, though they 
are consistent within the large uncertainty from pDIS. 
For GS-C and DSSV, which have a node in ^g{x) near 
the center of the sampled x region, a cancelation be- 
tween the positive and negative contribution in the wide 
X distribution in each px bin leads to a small value of 
^^([0.02,0.3] ^j-^j thus small All- As these two parame- 
terizations have significantly different AG'°'^' values, it 
is clear that we are sensitive only to AG^^ ^^'*^ '^! . 

For each pDIS fit, by varying AG^"'^' at the input 
scale while fixing the quark distributions and the shape 
of Ag(x) to the best fit values, a range of All curves 
were generated. Figure [3j3 shows the resulting x^ pro- 
files. While this approach is different from that in Fig.[2l3, 
resulting in a different profile shape for GRSV, the 
Ax^=l and 9 constraints are quite consistent. The 

^(^[0.02,0.3] ^j^i^gg ^jjg 

X^ minimum for each param- 
eterization are between 0.1 and 0.2, and are listed in 
Table [T| At Ax^=9, the profiles are consistent for all 
parameterizations, independent of shape, indicating that 
our data are primarily sensitive to the size of AG^'' ''^'^'^^ . 
The cross section for tt" production has been pre- 
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FIG. 4: (color online) profile as a function of AG^°-°l^-^^ 
when the theoretical scale is set to fj, — pt, Pt/2 and 2pT. 

sented Q and compared with NLO pQCD expectations 
with the theoretical scales (factorization, fragmentation, 
and renormalization) in the calculation all set equal to 
fi = kpT^ with k — 1. The calculation agreed with our 
results within the sizable theoretical uncertainties in the 
choice of scale, which were estimated by varying k up 
and down by a factor of two. As we rely on NLO pQCD 
to extract AG'[° °^'°''^1 from our measured ^lL' must 
consider the effect of this uncertainty. Figure |4] shows the 
change in the AGg^gy '^' constraint when varying k in 
the ^LL calculation in the GRSV parameterization. The 
theoretical scale uncertainty for the constraint on posi- 
tive values of AGq^Iy "^' is similar to that for varying 
the parameterization, while large uncertainty arises for 
negative values. 

We have presented results for A^^^ from Run-6, which, 
combined with Run-5 results, at /i^ = 4 GeV^ give 

^<^GRSV^' = 0.2 ±0.1(stat)±0.1(sys) 

t^;!^ (shape) ±0.1 (scale). (2) 

Using four parameterizations of A 17(2;), we find a shape 
independent constraint of -0.7 < AGI^ ^^^^-^I < 0.5 at 
Ax^=9 (~3fT). The theoretical scale induced uncertainty 
is small for positive values of AG^j^gy '^^ , but is sizable for 
negative values. Future measurements will be required 
to measure ^g{x) for x < 0.02 where large uncertainty 
remains and which may still contribute a significant 
amount of the proton spin. The quark spin contribution 
was well constrained by pDIS, and our result begins to 
significantly constrain the gluon spin contribution as well. 
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